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Abstract 



We show that the infinite imaginary anomalous dimensions of certain heavy quark- 
antiquark currents in the HQET arise due to an inappropriate commutation of two limits, 
namely coinciding velocities and infinite cut-off. This commutation is not apparent when 
dimensional regular ization is used, but it can be made manifest in a cut-off regular ization. 

We argue that operators containing heavy quark and antiquark fields with the same 
velocity in the HQET must not be matched to perturbative QCD but to Coulomb-type 
bound state contributions. We show how to do so at one loop level. 
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1. Introduction 

The so called Heavy Quark Effective Theory (HQET) has become a standard tool 
to study hadrons composed of a single heavy quark [1] (see [2] for reviews). For these 
hadrons there are only two relevant scales, tuq the mass of the heavy quark and Kqcd- If 
rriQ is large enough one can use perturbative QCD to calculate contributions coming from 
momenta p, p » \i » Aqcd where \x is a cut-off such that txiq » \i. For momenta 
p « \x one may not carry out reliable perturbative calculations but one can use HQET, 
the symmetries of which relate different form factors and matrix elements. Since there is 
no relevant physical scale between rriQ and Aqcd, one can match HQET to QCD at the 
scale rriQ order by order in perturbation theory. The Wilson coefficients thus obtained 
may be run to lower scales by using standard renormalization group techniques. Roughly 
speaking, for the above mentioned hadrons a current in full QCD can be written as a 
current in the HQET multiplied by a Wilson coefficient which depends on «(mg), ck(aO 
and an anomalous dimension which is computable perturbatively. 

It was initially thought that the HQET had no applications to hadrons containing 
more than one heavy quark, in particular to heavy quarkonium. The reason is that those 
hadrons can be understood as weak coupling Coulomb Type Bound States (CTBS), which 
cannot be obtained from the static approximation. Unlike in the heavy-light systems, the 
short distance effects are dominant in heavy quarkonium. Long distance (non-perturbative) 
effects have been traditionally simulated by the inclusion of a smoothly varying background 
field which eventually gives rise to contributions depending on the gluon condensate [3,4] 
(see [5] for a recent work). Although this method encodes important features of the 
non-trivial QCD vacuum, it is not obvious that it exhausts all possible non-perturbative 
contributions*. 

Surprisingly enough, the HQET seems to know about quarkonium. On the one hand, 
the infinite imaginary anomalous dimensions found in operators with a heavy quark and 
a heavy antiquark when the velocities coincide have been related to the static interquark 
potential [7]. On the other hand, quarkonia states appear to be related to Goldstone modes 
arising from the spontaneous breaking of U(ANhf) — > U(2Nhf)®U(2Nhf) in a HQET con- 
taining quarks and antiquarks with the same velocity [8,9]. Recently, it has been shown 
that HQET techniques can indeed be used to parametrize certain non-perturbative contri- 
butions in heavy quarkonium [9] . These contributions correspond to quantum fluctuations 
that put quark and anti-quark in quarkonium almost on shell when the relative three mo- 
mentum is much smaller than the inverse Bohr radius. In this kinematical situation the 
quark and the antiquark move almost with the same three velocily, namely the quarko- 
nium three velocity. Non-perturbative contributions in this region are important and can 
be described using HQET. However, in this situation the HQET must not be matched to 
QCD order by order in perturbation theory. One should keep in mind that in addition to 
rriQ and Aqqd a relevant scale, namely the Bohr radius a n ~ n/(arriQ), laying between 

* That would be the case if some non-perturbative contributions could not be expressed 
in terms of local condensates. This situation has been recently discussed in [6]. 
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rriQ and Aqcd enters the problem. The HQET here must be matched to CTBS contri- 
butions, which arise from an infinite summation of one gluon exchange diagrams of QCD, 
the so-called ladder diagrams. The matching above has been carried out at tree level in 
[9]. 

The main goal of this note is to explain how such a matching can be carried out at 
one loop order. The fact that HQET must be matched to CTBS contributions and not 
to perturbative QCD is related to the infinite imaginary anomalous dimensions found in 
[7,10]. We shall also see that those are an artifact due to an inappropriate commutation 
of two limits, namely coinciding velocities and infinite cut-off. 

In order to avoid confusion, we would like to properly locate our work with respect to 
the standard approach by effective theory techniques to heavy quarkonium, namely the so 
called Non-Relativistic QCD (NRQCD) [11]. NRQCD is an effective theory for momenta 
of the order of the inverse Bohr radius or smaller, in contradistintion to HQET which is 
an effective theory for momenta of the order of Aqcd- NRQCD also admits an expansion 
in 1/rriQ but momentum counts as inverse Bohr radius rriQa whereas energy as rriQa 2 . 
There is no flavour symmetry although spin symmetry remains. Contributions coming from 
momentum larger than the inverse Bohr radius can be calculated perturbatively in a(mo). 
They give rise to Wilson coefficients when matching NRQCD to QCD. Momenta of the 
order of the inverse Bohr radius give rise to weak coupling CTBS. However, there are also 
momenta of the order of Aqcd which bring in non-perturbative contributions. One may 
try to evaluate both perturbative and non-perturbative contributions by lattice simulations 
starting from NRQCD [12]. We believe, however, that one can proceed analytically one 
step further in the separation of perturbative and non-perturbative contributions, at least 
for very heavy quarks. We have been able to describe some of these non-perturbative 
effects by using HQET. Here, we would like to present in somewhat more detail how the 
matching of HQET to CTBS contributions can be carried out. It should be clear that this 
is quite different from describing quarkonia states starting from the HQET (static limit) 
which is known to be a hopeless task. 

We organize the paper as follows. Section 2 contains a technical discussion on the 
integral from which the infinite imaginary anomalous dimensions arise. In Section 3 we 
sketch how the HQET can be matched to CTBS. The result depends crucially on the 
correct evaluation of integrals analogous to the one analysed in Section 2. In Section 4 we 
carry out the matching at one loop within a proper setting. Section 5 is devoted to the 
conclusions. 

2. Dimensional regularization versus three-momentum cut-off regulariza- 

tion 

Consider the following integral, which arises in the calculation of anomalous dimen- 
sions of heavy-heavy currents in the HQET. 

I ± (v.v', v.k, v'.k') := -i J ( ^ 4 ^ + * fc) - ^ ±v , {q I ¥) +ie / +ie . (2-1) 
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v° , v ' >0. 

I± above is ultraviolet divergent. If we use dimensional regularization, i.e. I± — > I±, we 
obtain 

I+(v.v', v.k, v'.k') = - — ~p 4~~2#( u '-' y ) + finite terms, (2.2) 

11 7T 

I® (v.v', v.k, v'.k') = — -—r(g(v'.v) — i—) + finite terms v^v', (2.3) 

4 — D 4tt £ 

4tt 2L 4-7> v.{k - k') + ze \v.k\ \v.k'\ 

+ v.k9(v.k)iir - v.k'e(-v.k')iir)}, (2.4) 

where 

£ := ^{V.vY - 1 , ^) : =^ ln ^|- ( 2 - 5 ) 

The anomalous dimensions of certain currents may be obtained from the coefficients of 
the divergent terms in (2.2)-(2.4). We have calculated (2.2)-(2.4) directly in Minkowski 
space by carrying out first the integral over the zero momentum component. No Feynman 
parametrization or similar tricks which may lead to ill-defined expressions have been used. 
Keeping track of the ie is crucial to avoid ambiguities. Notice that 

lim I°(v.v', v.k, v'.k') ^ 7^(1, v.k, v.k') , (2.6) 

v' — >v 

which is going to be important in the following. The lhs becomes infinite imaginary when 
v — > v' (i. e. £ — * and g(v'.v) — > g(l) = 1). However, the rhs gives a finite real anomalous 
dimension. The original way of obtaining I®(v.v' , v.k, v'.k') [7,10] was through the analytic 
continuation of I+(v.v', v.k, v'.k') to negative v ' , more precisely v.v' — > —v.v' — ie , which 
coincides with the result (2.3). This procedure requires v ^ v' and hence the discontinuous 
limit (2.6) cannot be uncovered *. 

It is difficult to believe that physics is discontinuous when v' is close to v so the 
question arises which expression (2.3) or (2.4), if any of them, is relevant in that situation. 
In order to gain a better understanding of the non-commuting limit (2.6), let us consider a 
spherically symmetric three-momentum cut-off regularization, i.e. I± — > I± where f d 4 q — > 
/ d(v.q) J M d(e\q) (e\v = 0) . We obtain 

I+(v'.v, v.k, v'.k') = —-^-^g(v'.v) ln^ + finite terms (2.7) 



* Note that the analytic continuation v.v' — > —v.v' — ie transforms one HQET (two 
heavy quark fields with velocities v and v' respectively) into a different HQET (a heavy 
quark fields and a heavy antiquark field with velocities v and v' respectively) . This is quite 
different from the usual analytic continuation which relates the s and t channel within a 
given relativistic QFT. In fact, it only makes sense if we regard both HQETs as emanating 
from a single QFT, namely QCD, where crossing holds. 



4 



and hence the anomalous dimension agrees with (2.2). The ultraviolet behaviour of 
turs out to be qualitatively different in the following two situations 




a) 



This can be seen from the expression 



, , v'.k' 

» \v.k 

v.v 

v'.k' 



v.v 
v.v' 



<< \v.k - 



v.v' 



(2.8) 



lf^(v'.v, v.k, v'.k') = 



1 1 



x In 



— + it 



dq< 



+ 



q + v.k + ie q — v.k — 



1 



M. + a + vL*L 
In v ' v ' v ' v ' 



ie 



— ie 



(2.9) 



A + v .k - ^ + ie q + v.k + ie _ J2_ + q + vLhL _ ie 

v.v' v.v' v.v' ^ v.v 

Recall that always v'.k', v.k « fx. In the case i) we obtain 



1 



i 7T 

I^Ou'.v, v.k, v'.k') = lu.fi — -(g(v'.v) — i—) + finite terms, 

Atx z £ 



(2.10) 



which leads to the same result as (2.3) for the anomalous dimension. In the case ii) we 
obtain 



It{v'.v, v.k, v'.k') 



1 



— -(2/x + v.k In - v'.k' In — 

' '■' '< - ' \v.k\ \v'.k' 1 



4it 2 v.k — v'.k' + it 
+ v.ke(v.k)in - v' .k'9(-v' .k')in) 



(2.11) 



When v — * v', £ — > we are always in the case ii) for any finite cut-off fx. From the 
logarithmically divergent terms in (2.11) we identify the same anomalous dimension as 
in (2.4), but there is an additional linearly divergent term. This term cannot appear in 
(2.4) due to the dimensional regularization prescription of putting scale invariant integrals 
to zero. From (2.8) and (2.11) it is apparent that the relevant result for 7_ when the 
velocities of the quark and antiquark are close has to do with (2.4) and not with the limit 
v' — > v of (2.3). The point is that a physical scale fi must definitely be introduced to 
separate the momenta from which we can start using HQET. With v ^ v' there are two 
scales, namely fi and fi(v.v' — 1). Both scales are large for fi large as far as v ^ v' . In 
that case dimensional regularisation and the three-momentum cut-off regularisation are 
equaly convenient and they give the same result for the anomalous dimension (which may 
be imaginary). However, when v' — > v the second scale goes from being large to vanishing. 
Since dimensional regularization is only able to provide a single scale (the subtraction 
point), the only way it has to show us the brutal change in the second scale is by means 
of the discontinuous limit (2.6). The three-momentum cut-off is more convenient in this 
case. It allows us to examine the relevant limit in terms of physical scales (2.8) and to 
uncover the qualitative difference in the UV behaviour of 7_ for v' much different from v 
and for v' similar to v, that is the linear divergence in (2.11). 
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The linear divergent term in (2.11) has indeed far reaching consequences: it is non- 
local. The non-local cut-off dependence indicates that the HQET for quarks and antiquarks 
with the same velocity does not behave as a renormalizable quantum field theory, when 
we calculate Green functions involving currents with quark and antiquark fields In this 
situation the HQET would be of little use if we did not know the high momentum physics 
which lies beyond the cut-off. Fortunately, we do know the high momentum physics which 
lies beyond the cut-off: Coulomb-type bound states. It is not obvious, however, whether 
one can actually match these contributions in the HQET with contributions arising from 
CTBS. We devote the remaining sections to this goal. 

After having found a non-local linear divergence in (2.11) we may wonder at which 
extend the result that we obtain depends on the regularization that we use. For v = v' and 
taking the Coulomb gauge in (2.1) (i.e. l/(q 2 + ie) — > — l/(e\g) 2 ) the following statements 
are easy to check. For both U(N h f) and gauge invariant regulators of the kind 

1 1 D 2 

f(--k) (2-12) 



fv.D + ie fv.D + ie \i 2 

where f(x) — > strong enough when x — > oo and /(0) = 1, the non-local linear divergence 
and the pole in (2.11) remain. Regulators of the kind 

(2.13) 



iiv.D + ie ^ D _^_^n + ie 



which do not respect the U{ANhf) symmetry, give different branching point singularities 
depending on n instead of a simple pole. L r (4A^/ l /)-invariant 'regulators' of the kind 

(2.14) 



1>v.D + ie f(v. D _ fl (_^ n y j +ie 



do not regulate the integral. Neither a Pauli-Villars regularization for the heavy quark (an- 
tiquark) propagators does. Consequently the non-local linear divergence in (2.11) appears 
to follow from the use of a L r (4A^/ l /)-invariant regulator. Indeed, the three dimensional 
cut-off regularization respects the U(4Nhf) symmetry. We have checked that the Ward 
identities for the heavy quark-antiquark current are fulfilled at one loop level with no need 
to introduce any counterterm* . 

To summarize, for v' very different from v (in the precise sense (2.8)) we find the same 
imaginary anomalous dimensions as previously obtained [7,10]. For v' close to v (in the 



t This is not in contradiction with ref. [13] where it was found that HQET, including 
quarks and antiquarks with the same velocity, is renormalizable since the renormalization 
of currents was not considered. 

* The three dimensional cut-off regularization, however, breaks local SU(3) gauge in- 
variance (Slavnov- Taylor identities) which, nevertheless, can be restored by adding suitable 
(non-gauge invariant) local counterterms. We have checked this stament at one loop order. 
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precise sense (2.8)), we conclude that there is no infinite imaginary anomalous dimension 
as claimed in [7], but a non-local linear cut-off dependence. This non local dependence 
cancels against contributions coming from high energies (CTBS) once the matching has 
been performed properly, as we shall see in Sect. 4. 



3. Matching CTBS with HQET: getting the flavour. 

Consider the following Green function to one loop in the HQET. 

G r (ki,k 2 ) = J d A Xl d A x 2 e- ik ^- ik ^{b\T\h a Th\n)^ |0> , (3.1) 

a, b... = 1, ...Nhf are flavour indices (repeated flavour indices are not summed up). Colour 
indeces are not explicitly displayed in the colour singlet currents. Otherwise they will be 
denoted by i±, i 2 , ■■■= l....N c , N c being the number of colours. a±, a 2 are spinor indices. h+ 
annihilates quarks whereas h a _ creates antiquarks. h± have two independent components 
each. We combine them often into a four independent component field h a := h+ + h a _ 
(h± = P±h;p± = (1 ± ^)/2). Let us choose the Coulomb gauge {diB 1 = 0, x 1 := e l .x, 
B l := e\B, e\v = 0, i = 1, 2, 3). We obtain at one loop level (see Fig. 1) 

G r (k u k 2 ) =G° r (k u k 2 ) + k 2 ) + G^ih, k 2 ) + G^ih, k 2 ). (3.2) 

Gr{k!,k 2 ) =(p-Tp + ) a2ai 8 ili2 1 1 (3.3) 

v.k 2 + ie v.ki + te 

G^(k u k2) =^li(p-Tp+) aaai 6 ilia ( 1 ) 2 (3.4) 
TV v.k 2 + it v.ki + ie 

G?{ku k2) = ^_ Tp+)a2M _l_ ? _i_ . (3.5) 

gpfr, h) = - ^ 2 , {p _r p+)maA ^ _ _J_ _ _ J_ _ a€ . ( 3. 6 ) 

N 2 - 1 

The linear divergences in (3.4) and (3.5) are local and may be understood as mass renor- 
malizations. The linear divergence in (3.6) is non-local and appears from an integral 
analogous to the one discussed in the previous section. If we take the Lorentz-Feynman 
gauge additional logarithmic dependences in fx arise. The choice of the Coulomb gauge 
is not only a matter of convenience. As we have already mentioned we should match the 
HQET results to CTBS contributions. The latter are always calculated in the Coulomb 
gauge. Then, in order to match with the results obtained by the conventional treatment 
of CTBS, using the Coulomb gauge in the HQET is mandatory. 

In ref. [9] it was shown that, at leading order in //flab,n? &ab,n 

being the Bohr radius, 

G®(k[, k' 2 ) can be matched with 
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G T {puP2) := J ^i^e^^+^^^OlTlQ^g^O)^ 1 ^!)^?^)} |0) . (3.7) 

Pi = -m b v H £/ a b,ni> - «i , = -m a v H E abtn v - k 2 , 

m ab m ab 

rriab :=m a + m b , k[ , k 2 — > , 

under the assumption that CTBS dominate the above Green function. The precise relation 
reads 

Gr(pi,P2) = C r G r (k[,k' 2 ) , C r = *: M (0)* aM (0). (3.8) 

E a b,m ^ab,n(x) and ^ a b,n(k) are the energy, the coordinate space wave function and the 
momentum space wave function respectively of a Coulomb-type state with principal quan- 
tum number n. We call (3.8) tree level matching, meaning that the CTBS contribution 
matches with a tree level diagram of the HQET. 

Let us concentrate on G]?(k 1 , k 2 ). Anologously to (3.8), we would like to find a 
contribution coming from CTBS which matches with this contribution in the HQET. In 
order to do so, one way or another, we have to single out a gluon propagator from the 
infinite set in the ladder diagrams. At the risk of having some double counting, which we 
shall overcome in next section, we may proceed as follows. We start as if we computed 
(3.7) at one loop in QCD (or NRQCD), namely 

Gf{PuP2) = J d 4 x 1 d 4 x 2 e tp ^ +tp2X2 (0\T{Q a TQ b (0) 

x J d 4 zd 4 z'i 2 {Q h 1 ^Q h ){z){Q a lv B v Q a ){z') ^ 

x<5£WQS a W}|o>. 

Pi and^2 are as in (3.7). At this stage we introduce resolutions of the identity approximated 
by CTBS between the currents and the interaction terms. Notice that since we use the 
Coulomb gauge the two interaction currents become a single time four quark operator. We 
then have (see Fig. 2) 

G 1 J ?(p 1 ,p 2 ) = j d 4 x 1 d 4 x 2 e tpiXl+ip2X2 J d 4 zd 4 z , i 2 e(-max(z°,z , °)) 

xe(rmn(z ,z'°)-max(xlx° 2 ))J2 / jSf^o j ( 2 ^5» (0\Q a rQ b (0)\n) 



n,m 

\brnr„, f^b 



x (n|T{(g 6 T^ M g 6 )(^)(g a TVQ a )(^0}k)(^inQa\H^i)^ 2 (^2)}|o) 

^ d 3 q i 



xg 2 



We always calculate in the non-relativistic limit a a b )TJ m a , a a fc )n m& >> 1. In this limit 
7a* ~ v fj,i which we have already assumed when we wrote down the static component of 
the gluon propagator only. If we constrain the gluon momentum to be smaller than //, we 
obtain at leading order in fxa a b, n (l^a a b,n << 1) 

Gi?( Pl ,P2) = - *: t ,„(0)*„6,„(0)-^2 /J (p_r p+ )„ aai 

i i i ( 3 -") 

x Si 



v.k[ + ie v.k' 2 + ie v.(k[ + k' 2 ) + ie 



which matches G}?(k' 1: k 2 ) in (3.6) with the same Wilson coefficient Cr as in (3.8). This 
is very illustrative. It tells us that the HQET at one loop reproduces the exchange of a 
singled out soft gluon with momentum lower than \x in a CTBS. If we add to (3.10) the 
momenta larger than \i for the singled out gluon the result, of course, is fx independent. 
Thus we have identified the large momentum contributions which cancel out the worrysome 
non-local linear divergence in G]?(ki 7 fo). 



Although the exposition above contains the essence of the matching between CTBS 
and HQET, the procedure has to be refined. There are indeed two outstanding problems. 
First of all, Gj, 1 (ki, k 2 ) and G r 2 (/ci, k 2 ) cannot be matched with (3.7). These contributions 
correspond to one particle reducible diagrams and amount to mass renormalizations of the 
heavy quark and the heavy antiquark respectively and independently. Contributions like 
these do not arise from a calculation of the kind (3.7)-(3.11), that is, from the diagrams 
Gp 1 (pi,P2) and G]?(pi,p 2 ) in Fig. 2. The reason is that inside quarkonia quark and 
antiquark are always communicating through the exchange of gluons and hence one cannot 
get independent contributions for either. We can get ride of this missmatch by restricting 
our analysis to local operators of the heavy quark and antiquark fields, as we shall do in the 
next section. This is not an important restriction since local operators are the objects which 
actually arise in the calculations of physical observables. Second, and more important, 
there may be a problem of double counting. On the one hand, we sum up an infinite set of 
gluon propagators of arbitrary momentum in ladders to obtain the Coulomb potential. On 
the other hand, we single out a gluon propagator, which we force to have small momentum, 
and treat it perturbatively. (Is this gluon propagator not already counted in the infinite 
set?). We show how to avoid this problem in the next section. 

4. Matching CTBS with HQET: a proper setting 

Consider the Green function analogous to (3.1) but with quark and antiquark fields 
at the same point in a current. 

Grr'(p) := J d 4 xe ipx (0\T {Q a TQ b (0)Q b T'Q a (x)} |0> , (4.1) 

p = -m ab , n v -k , k -> . 
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We introduce two cut-offs /j/ and fx, both of them much smaller than the inverse Bohr 
radius. The first cut-off \i' separates two different kinematical situations, small (high) 
relative three momentum, in the heavy quarks currents. We denote them as on-shell (off- 
shell) currents. Our goal is reproducing (4.1) in the HQET when one of the currents is 
almost on-shell. 

G° r n r(p) ■■= f d A xe^(0\T{(Q a TQ b ) off (0)(Q b T'Q a ) on (x)}\0), (4.2) 

The on-shellness will be implemented by cutting off the relative momenta larger than // in 
the wave function. The second cut-off fi separates the gluon momenta between large and 
small. We assume that the relevant contributions from large gluon momenta are the infinite 
sums of one gluon exchange diagrams (ladder diagrams) which give rise to the CTBS. For 
small gluon momenta we use perturbation theory (no infinite summation) t. Again we 
work at leading order in \ia a b,n and in the non-relativistic limit Tn a a a i, jril rnb&ab,n 

» 1. 

Notice that double counting is explicitely avoided by the introduction of the cut-off. This 
is a proper framework in which the matching between CTBS and HQET can be carried 
out. A given order of perturbation theory for small momentum gluons calculated this way 
corresponds to the same order of perturbation theory in the HQET. We carry out the 
matching at one loop below. 

At zeroth order, i.e. no soft gluon propagator, we have 



G&(p) = -,NM P -T P+ r,K^ K , nm £) - - J ab M - u , (4.3) 

AE ah ^ n (p) := E ab , n ([i) - E abjrL . 

Because we only include momenta larger than fx in the gluon propagator giving rise to the 
Coulomb potential, both the wave function and the bound state energy depend on /i. 
More precisely the Coulomb potential is modified by 



AV(p,i?) = C f g 2 9(v - \p- f?\) 1 . (4.4) 



Since [ia ab ^ n « 1, AV can be treated as a perturbation. We may evaluate ^>^ b n and 
AE a (, in (/i) at any desired order in quantum mechanical perturbation theory. 

In the HQET we have 

Grr'(fc) := J d 4 xe~ lkx (0\T {h a _Th b + (0)h b + T'h a _(x)} |0> . (4.5) 



t For rriQ — > oo we can always assume that perturbation theory is applicable at scales 
even lower than rriQa since this quantity also goes to infinity and a(niQa) may still be 
small. This is not so for the actual values of the nib and m c . Nonetheless once we find the 
structure of the matching perturbatively we believe we may safely assume that the same 
structure holds when one takes into account non-perturbative contributions as well. 
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At one loop 

Grr'(fc) = G rr °(k) + Grr, 11 ^) + G rr > 12 (k) + G rr > 13 (k) . (4.6) 

where Grr ,0 (k) stands for the tree level contribution. It is easy to see that (4.3) matches 
Grr" (&) with the following Wilson coefficient Cp(/f) and modification of the HQET la- 
grangian AL 

Cv{n) = *:i(0)^ 6) n(0) 

AL = l -AE ab Mh a h a + ±AE ab , n (,u)h b h b . (4.7) 

That is 

Grr'(p) = C r (iJ.)Grr'(k). (4.8) 



With the modifications above the HQET lagrangian gives rise to the following one 
loop contributions (see Fig. 3) 

- 2Grr' 11 (£0 = -2G r r' 12 (£0 = G r r' 13 (£0 = 

- iiWMwrO(£)( ..,. * ,.,.-. )'SgV <4 ' 9) 



67T 2 v.k + AE ah ^ n (p) + ie' 



n 



The contributions of the full theory (4.2) at one loop can be written as (see Fig. 4) 
G°MP) =G^\p) + G^\ P ) + G r n r/ 12 (p) + G r n r , 13 (p) (4.10) 
Gft, n (p) := J d 4 xe tpx (0\T{Q a TQ b (0) J d 4 zd 4 z'- 



,i 2 



2 



x (Q a 1 ^Q a )(z)(Q a 1 ^Q a )(z')Q b T'Q a (x)}\0) . (4.11) 
G° T n T , 12 {p) :=J d 4 xe ipx {0\T{Q a TQ b (0) J d 4 zd 4 z' 1 - 

x (Q b ^Q b )(z)(Q b ^Q b )(z')Q b T'Q a (x)}\0) . (4.12) 
G™, 13 (p) := J d 4 xe ipx (0\T{Q a TQ b (0) J d 4 zd 4 z'i 2 

x (Q b ^Q b )(z)(Q a ^Q a )(z')Q b T'Q a (x)}\0) . (4.13) 

In order to account for momenta larger than \i we insert resolutions of the identity between 
the currents, which we approximate by CTBS, as we did in (3.10). Recall that the gluon 
propagator must be cut-off at \i « l/a a b, n and that almost on-shell contribution is 
implemented by cutting off the relative momenta larger than fx' « l/a a & ;n . We obtain 



2Gft, n (p) = -2G° r n r , 12 (p) = G° r n r , 13 (p) = 

/3 



1 .oCfC* (4-14) 
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It is remarkable that each of the contributions above matches the corresponding one 
loop contribution in (4.9) with the same modification in the HQET lagrangian and the 
same Wilson coefficient (4.7). That is (4.8) holds at one loop as well. 

Although the matching stated as above makes perfect mathematical sense as it stands, 
it still requires some adaptation to the precise physical situation we are concerned with. 
Namely, in a CTBS we wish to separate large and small gluon momentum contributions. 
A CTBS is obtained through an infinite sum of one gluon exchange diagrams only. Con- 
sequently in the perturbative evaluation of Grr'ip) we should only keep diagrams which 
correspond to one gluon exchange. This means, for instance, that at one loop (4.11) and 
(4.12) must be dropped, since these contributions correspond to quark and antiquark self- 
energies. Once this is done, it is easy to see that Grv{p) becomes \i independent if we 
calculate n and AE ab)n ([i) at the same order in quantum mechanical perturbation 
theory as we have calculated the small momentum contributions in perturbation theory. 

For instance, since we used zeroth order in perturbation theory in (4.3) we have to 
calculate ^>^ b n and AE ab ^ n (fi) at zeroth order in quantum mechanical perturbation theory. 
That is 

K b ,n = *a b ,n , AE ab ,M = , (4.15) 

which substituted in (4.8) gives the result found in [9]. At one loop order we must calculate 
^ab n an d AE abjn (/i) at first order in quantum mechanical perturbation theory. We obtain 

C ex 

Kb,n = *ab,n , AE ab ,M = f —2». (4.16) 

If we next consistenly expand G rr ,°(p) in a\i at first order we see that the second term in 
the expansion exactly cancels Gpp, 13 (p). Since we dropped (4.11) and (4.12) we must also 
drop their counterparts in the HQET in order to keep (4.8) true. Fortunately, this can 
be safely implemented by the addition of the following mass counterterms in the HQET 
lagrangian t 

All = -^yTh a h a - ^iih h h b . (4.17) 

It is again remarkable that the mass counterterms above exactly cancel the effective 
masses —AE ab ^ n (ii) /2 in (4.7) at this order. As the total outcome we have that at one 
loop the Wilson coefficient is the same as at zeroth order and that the HQET lagrangian 
must not be modified either. It is not clear to us whether the fact that neither the Wilson 

t (4.17) has a clear physical meaning. It arises because when integrating large momenta 
we do it from an infrared cut-off \x onwards. After doing so, the pole mass depends on 
\x. However, we do not actually work with a \x dependent mass but with the pole mass at 
[i = 0. The ^-dependent pole mass and the pole mass we work with differ by an amount 
5m ~ \i. This difference is the origin of (4.17). Notice that this is much in the same way 
as we use -Z? a 6,n instead of E a b^ n {\i) and hence the effective mass term (4.7) arises. 
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coefficient nor the HQET lagrangian receive corrections at one loop remains true at higher 
orders in perturbation theory. 

The maching we have carried out is exact when tuq — > oo. The philosophy is analogous 
to the matching carried out in heavy-light systems. However, the matching here is done 
in perturbation theory for small momentum gluons only. Large momentum gluons are 
summed up in the ladder diagrams in order to account for the CTBS. Then the Wilson 
coefficients have to do with the wave functions instead of with anomalous dimensions. 

Once we have established in the proper framework above the matching between HQET 
and the CTBS, we can safely use the symmetries of the former in order to get information 
in a region where perturbation theory is not applicable, as it was initially done in [9] . 

5. Conclusions 



We have pointet out that there are no infinite imaginary anomalous dimensions in the 
HQET. When the quark and the antiquark fields have the same velocity the HQET must 
not be matched to perturbative QCD but to CTBS contributions. We have shown how to 
do it at one loop level. 
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Caption 1. Diagrams contributing to Gr at one loop in the HQET theory. The 
heavy quark field carries residual momentum k\ and the heavy antiquark field residual 
momentum h^. 

Caption 2. Diagrams contributing to Gr at one loop in the full theory. The quark 
field carries momentum p\ and the antiquark field momentum pi- The double line means 
almost on-shell, the wavy line means soft gluon and the square means summing up all the 
ladder diagrams. 

Caption 3. Diagrams contributing to Grr' at one loop in the HQET theory. The 
heavy quark anti-quark current carries momentum k. 

Caption 4. Diagrams contributing to G r p, at one loop in the full theory. The current 
carries momentum p. 
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